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Abstract 

We calculate the temperature-dependent charge carrier transport of bilayer graphene 
(BLG) impacted by Coulomb impurity scattering within the random phase approxima- 
tion. We find the polarizability is equal to the density of states at zero momentum 
transfer and is enhanced by a factor log 4 at large momentum transfer for arbitrary 
temperature. The sharp cusp of static polarizability at q = 2kp, due to the strong 
backward scattering, would be smooth by the increasing temperatures. We also obtain 
the asymptotic behaviors of conductivity of BLG at low and high temperature, and find 
it turns from a two dimensional electron gas (2DEG) like linear temperature metallic be- 
havior to a single layer graphene (SLG) like quadratic temperature insulating behavior 
as the temperature increases. 

PACS number(s): 81.05.Uw; 72.10.-d, 72.15.Lh, 72.20.Dp 

1 Introduction 

Since graphene, a two-dimensional single layer of graphite, is fabricated [1] which has at- 
tracted much attention from both experimental and theoretical physicists. One important 
experimental puzzle is that there is a so-called "minimum conductivity" at the charge 
neutral (Dirac) point. Several early theoretical work [2J calculated a universal T = min- 
imum conductivity a m i n = 4e 2 /7r/i at the Dirac point in disorder-free graphene, but the 
experiments show that the conductivity has a non-universal sample-dependent minimum 
conductivity plateau (~ 4e 2 //i — 20e 2 /h) around the Dirac point. By using a random-phase 
approximation (RPA)-Boltzmann formalism, this puzzle has been theoretically explained 
as the result of carrier density fluctuations generated by charged impurities in the substrate 
[3j HI [5] . Therefore, the screened Coulomb scattering plays an important role in understand- 
ing the transport properties of graphene, and several corresponding theoretical researches 
P El El [9] have been made. 

While single layer graphene (SLG) has been widely studied from both experimental and 
theoretical sides, bilayer graphene (BLG), which is formed by stacking two SLG in Bernal 
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stacking, as an other significant carbon material, is attracting more and more attentions \10\ 
QTJ [12] due to its unusual electronic structure. BLG has a quadratic energy dispersion [10] 
similar to the regular two dimensional electron gas (2DEG) but its effective Hamiltonian is 
chiral without bandgap similar to the SLG. Although it is reported [13J that a widely tunable 
bandgap has already successfully realized in BLG by using a dual-gate bilayer graphene 
field-effect transistor and infrared microspectroscopy, here we still ignore the bandgap in 
BLG dispersion and keep the transport properties of such BLG with tunable bandgap as a 
question studying in other paper. 

Because of the role of screened Coulomb scattering by charged impurities in understand- 
ing the transport properties of SLG, it is significative to investigate the affection of screened 
Coulomb scattering by charged impurities in BLG. The screening function of BLG at zero 
temperature [14J and SLG at finite temperature [15] have already been analytic investigated, 
however, the BLG analytic form at finite temperature and the corresponding temperature- 
dependent behaviors of transport, which are the issues we will discuss in this article, has 
not yet been provided. Although it is argued in Ref. |16j that any strong screening-induced 
temperature dependence should not been anticipated in BLG resistivity and the relatively 
strong collisional broadening effects would suppress the small screening-induced temper- 
ature dependence due to the small mobilities of current bilayer graphene samples, it is 
still significative to investigate such screening-induced temperature-dependent behavior for 
comparing the affections of screened Coulomb scattering by charged impurities on transport 
properties in BLG, SLG and 2DEG and representing how the BLG behaves as the crossover 
from SLG to 2DEG. 

This article is organized as the following. In Section 2, we present the Boltzmann 
transport theory to calculate temperature-dependent bilayer graphene conductivity. In 
Section 3, the temperature-dependent screening function is investigated. In Section 4, we 
present the asymptotic behavior of conductivity at low and high temperature, and numerical 
results obtained. The conclusion is given in Section 5. 

2 Conductivity in Boltzmann Theory 

The BLG Hamiltonian has an excellent approximate form in the low energy regime which 
can be written as [10] (we set Ti = 1 in this article) 



where m = r y\/{2v F ) 0.033m e , 71 0.39e^ is the interlayer coupling, and vp 10 6 m/s 
is the SLG Fermi velocity. The corresponding eigenstates of Eq.(pQ) are written as 
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where L 2 is the area of the system, s = +1 and —1 denote the conduction and valence 
bands, respectively, and 9^ = arctan^ / k x ) is the polar angle of the momentum k. The 
corresponding energy is given by e g ^ = sk 2 /2m, and the BLG density of states (DOS) is 
iVo = mg/2-rr (g = g v g s = 4 is the total degeneracy.) which is a constant for all energies and 
wave vectors. 

When the electric field is small, the system is only slightly out of equilibrium. To the 
lowest order in the applied electric field E, the distribution function can be written as 
fsk = f( e sk) 9 k wnere /( e s fc) ^ s t ne equilibrium Fermi distribution function and g ? is the 
deviation proportional to E. The Boltzmann transport equation is given as 



dt 



dt dk 



-eE ■ v 



df_ 



d 2 k 



2 \9sk 9 S 



w 



sk,sk' ' 



where v ? = sk/m is the electron velocity, 
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n, is the number of impurities per unit area, and (V^ g) is the matrix element of scattering 
potential with an average over configuration of scatterers. For elastic impurity scattering, 
the interband processes (s ^ s 1 ) are forbidden. Under the relaxation-time approximation, 
we get 
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where the scattering time T(e s ^) is given by 
1 
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and 9^g is the scattering angle between k and k'. 
We know that the electrical current density is 

d 2 k 

3=9' 



Then we can get the electrical conductivity by using Eq.([6|) 

iV oe 2 



in 
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/(efc) is the Fermi distribution function /(e^) = {1 + exp[/?(efc — //)]} 1 where (3 = 1/ksT 
and /i is the finite-temperature chemical potential. At T = 0, /(e) = 9{ep — e) (where 

£f = //(T = 0)), then we get the conductivity formula a = ^ £L Nor(Ep) which has the 
same form as the usual conductivity formula. 

The matrix element of the scattering potential of randomly distributed screened charge 
impurity in BLG is given as 



(V, 



sk,sk' 
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where q =| k — k' |, 9 = and t>i(g) = 2ite 2 /{nq) is the Fourier transform of the potential 
of the charge impurity with a background dielectric constant k. The factor (1 + cos 29) /2 is 
derived from the sublattice symmetry of BLG, while this factor is replaced by (1 + cos#)/2 
for SLG. The finite-temperature RPA dielectric function can be written as e(q) = e(q, T) = 
1 + v c (q)U(q, T), where v c (q) is the Coulomb potential and Tl(q,T) is the irreducible finite- 
temperature polarization function. Then the scattering time for energy e& of BLG is written 
as 

1 f d 2 k' . Vi(q) 



™ l I (2^)2 1 7TT |2 6{ek " efc ' )(1 " cos ^ )(1 + cos2 ^' 



r(e fc ) 'J (2vr)2 1 £ {q) 
comparing to the scattering time of SLG 
1 _ f d 2 k' Vi(q) 



T ( e k) S LG J ( 27r ) 2 £ (?) 

and the scattering time of 2DEG 



p J2U ( \ 

nni J 1 7T) | 2 ^-^)(l-cos^)(l + cos0), (12) 



2 ^'. / ^\ V M\ 2 5{e k -c k ,){l-cos9). (13) 



T ( e k) 2 DEG " ^ ( 27r ) 2 £ (<?) 

We can find that formally the three formulas are almost the same except the angular 
factor which arises from the sublattice symmetry, (1 + cos2#)/2 for BLG, (1 + cos#)/2 
for SLG and 1 for 2DEG. Actually the dielectric function e(q) are also different for the 
three systems, which would lead to different scattering times in the three systems except 
for the angular factors. They all have the same factor (1 — cos 9) which weights the amount 
of scattering of the electron by the impurity and always exists in Boltzmann transport 
formulism. This factor (1 — cos#) favors large-angle scattering events, which are most 
important for the electrical resistivity of the regular 2D systems. However, in SLG the 
large-angle scattering, in particular the 2fcp backward scattering, is suppressed by the factor 
(1 + cos 9). In contrast to the regular 2D system, in SLG the dominate contribution to the 
scattering time comes from the /cp "right-angle" scattering (i.e. 9 = tt/2) but not the 
2/cf backward scattering. Different from SLG, the 2kp backward scattering of the BLG 
is restored and even enhanced by the factor (1 + cos 29) which arises from the sublattice 
symmetry of BLG. Because of the restoral of the 2/cp backward scattering, many theoretical 
approaches which fit the ordinary 2D systems can been used for the BLG. Due to the 
qualitative similar, in some regimes the temperature-dependent behavior of polarization 
function and the transport properties in BLG are more similar to the 2DEG than the SLG 
as we will show below. 

3 Temperature-Dependent Polarizability and Screening 

First let us consider temperature-dependent screening 

2vre 2 

e(g,T) = l + — n( 9 ,T), (14) 
nq 
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where H(q,T) is the BLG irreducible finite-temperature polar izability function, which is 
given by (calculated at T = in Ref. [13] for BLG) 



U(q,T) = L£—Lrt Faa ,{k,9), (15) 



here k' = k + q, = sk 2 /2m, and F ss r{k,k') = (1 + ss'cos26')/2 where 6* = 0£r>, / s g 
is the Fermi distribution function / r = [exp{/3(e s £ — /i)} + where /i = /i(T) is the 
finite-temperature chemical potential determined by the conservation of the total electron 
density as 



where Tp = ep/kB and 



It is easy to find that 



^=F (^)-F (-P fJ ,), (16) 



t n rlt 

Fn{x) = I ' - , (17) 

exp(i — x) + 1 



F (x) =log(l + e a; ), (18) 

substitute it into Eq.(16), then we obtain 

KT) = e F , (19) 

which means that the chemical potential of BLG is temperature-independent and very 
different from that of the SLG and the regular 2D systems. 
We rewrite the polar izability as 

U(q,T) = n mtra (q,T)+U inter {q,T), (20) 

H-intra an d Winter indicate the polarization due to intraband transition and interband tran- 
sition, respectively, which are given by 

IWg,T) = -^E ^:^ ,1 + C °f % ; (21) 
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and 



U mter (q,T) = -f^E J ~ sk (22) 

where e s £ = sk 2 /2m, k' = k + q, and 
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here <j) is an angle between k and q. After angular integration, we obtain 
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here An = mg/2ir is the BLG density of states, f(e) is the Fermi distribution function 
f(e) = [exp{(3(e — //)} + . Then we have the extrinsic BLG static polar izability at finite 
temperature as 



U(q,T) = A J^y^ + q^-k 



2 'k 2 -q 2 
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At high temperature (T > Tp), Eq. ([26|) can be written as 
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At low temperature (T <C Tp), Eq. (|26p can be written as 
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(for q > 2k F ), (29) 
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Figure 1: Temperature-dependent BLG polarizability (a) as a function of wave vector for 
different temperatures and (b) as a function of temperature for different wave vectors. Here 
-/Vq = mg/2ir. 



For q = 2k F , we have 



s^> K c _ ^ _ ^ (1/2) ( _l ) - _ ^ _ | K(3/2) ( - )- , (34) 



where C = \/5 — log[(l + \/5)/2], £(x) is Riemann's zeta function. From above, we can 
give the doped BLG polarizability at zero temperature which is same as the result firstly 
obtained in Ref. |14j . 

n(g '^ = 0) = g (q) - f (q)6(q - 2k F ). (35) 
iV 

The screened potential is 

tt( \ _ "(g) _ 27re 2 = 27re 2 
W e(g) K?[l + t; c II(<z)] *(g + g s ) 5 1 ' 

where q s (q,T) = qv c (q)U(q,T) = 27re 2 U(q, T)/k = qTF~R(q,T)/N with q TF = mge 2 / 'k 
being the Thomas-Fermi screening wave vector of BLG. It is interesting to see that in the 
q — > long wavelength limit, the q s (q, T) of BLG is a constant value for all temperatures, 

q s (q = 0,T) = q TF = Ar s k F , (37) 

which is remarkably different from that of the SLG (see Eq.(29) and (30) of Ref. |15]). 

The BLG finite-temperature polarizability H(q, T) as a function of wave vector for dif- 
ferent temperatures and as a function of temperature for different wave vectors are shown 
in Fig.l (a) and (b), respectively. One novel phenomenon is that at q = 0, the BLG po- 
larizability equals to a constant value for all temperatures, i.e., U(q = 0,T) = Nq. This 
is a qualitative difference between BLG and SLG (or 2DEG) polarizability function, while 
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Figure 2: Temperature-dependent SLG polar izability (a) as a function of wave vector for 
different temperatures and (b) as a function of temperature for different wave vectors. Here 
A> = gk F /(2-Kv F ). 




Figure 3: The 2DEG polar izability function ( function of wave vector for several 

different temperatures T/Tf = 0,0.2,0.4,0.6,0.8,1.0 (top to bottom), (b) as a function of 
temperature for different wave vectors. Here Nq = rag/ir. 
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the latter H{q = 0, T) as a function of temperature changes notably (as shown in Fig. 2 for 
SLG and in Fig. 3 for 2DEG). The reason is, at q = 0, the intraband transition polarization 
Hi ntra (0,T) = Nq while the interband transition polarization ILi nter (0, T) = 0, i.e., the in- 
terband transition is forbidden at zero momentum transfer for all temperatures. The other 
remarkable phenomenon is the BLG polarizabilty approaches a constant value Nq log 4 in 
the large wave vector regime, arising from the fact that the interband transition dominates 
over the intraband contribution in the large wave vector limit. Therefore the domina- 
tive contribution to the whole polarizability has a crossover from intraband transition to 
interband transition for all temperatures. The weak temperature-dependent behavior of 
polarizations in the small and large wave vector regimes is a distinctive electronic property 
of BLG. Different from that of BLG, as shown in Fig. 2(a), the polarizability of SLG in- 
creases monotonically for large q stemming from the domination of exciting electrons from 
the valence band to the conduction band while the polarizability of 2DEG in large q limit 
decreases as 1/q 2 (as shown in Fig. 3 (a)). 

In contrast to the SLG and the 2DEG, the polarizability of BLG as a function of wave 
vector shows a nonmonotonicity, i.e., the polarizability at T = is monotone increasing 
with q in the regime [0, 2/cp] and monotone decreasing in the regime larger than 2kp. In 
ordinary screened Coulomb scattering, the most dominant scattering happens at q = 2kp, 
which gives rise to the famous Friedel oscillations. Due to its sublattice symmetry, the 2kp 
backward scattering of SLG is suppressed and therefore there is no singular behavior happen 
at q = 2kp. However, in BLG the 2kp backscattering is restored and even enhanced because 
of its chirality, which leads to a sharp cusp and a discontinuous derivation of polarizability at 
T = 0. We find that the temperature dependence of BLG polarizability at q = 2kp is similar 
to that of 2DEG polarizability, both are much stronger than that of SLG polarizability. Due 
to the strong temperature dependence of the polarizability function at q = 2kp, as shown 
below, the BLG would have a anomalously strong temperature-dependent resistivity for 
T -C Tp which is similar to that of the regular 2D systems [T7j. We also find the strong 
thermal suppression of the singular behavior of polarizability at q = 2k p, which is similar 
to that of the 2DEG. 

Showing the difference between BLG and other 2D systems, we provide the polarizability 
functions of SLG and regular 2DEG in the regimes of q = and q = 2kp in the low (T <C Tp) 
and high (T 3> Tp) temperature limits. For T <C Tp, 
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{for SLG) 



(39) 



{for 2DEG) 



here -Do = 9Ef/2t:v f and Nq = gm/2'K are the density of states of SLG and regular 2DEG 
at Fermi level, respectively. Comparing to the corresponding screening formula for BLG, 



n(g = 0,T) 
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U(g = 2k F ,T) « No [C - ^(l - V2)a\){^r) 1/2 ]. 



(41) 



We can find at q = 2k F the zero-temperature value of polarizability (normalized to the den- 
sity of states at Fermi level) of BLG is different from that of 2DEG, but their temperature- 
dependent parts are both the same, which represents the similarity of BLG and 2DEG. 
For T > Tp, 



F 



n(g,T) 



1 



6fef T 



(for SLG) 
(for 2DEG) 



(42) 



The corresponding high-temperature screening formula for BLG is given by 

q 2 T F ' 



U(q,T)^N 



1 + 



6k 2 F T 



(43) 



In high temperature limit, the polarizability of BLG approaches a constant value (i.e. Nq), 
which is very different from that of SLG, where the static polarizability increases linearly 
with T, and the regular 2DEG, where the polarizability falls as 1/T. The BLG shows an 
intermediate behavior between the SLG and the regular 2DEG. 



4 Conductivity Results 

4.1 Analytic Asymptotic Results 

In this section, we study analytically the static conductivity of the BLG in low and high 
temperatures limit. Firstly let us consider the temperature dependence of conductivity in 
the low temperature limit (T <C T F ). Using Eq. (fTTI) . the scattering time t(e f ,T) at the 
Fermi level ep in the Born approximation is given as 

n t f 2k ^ r q 2 [l - 2(q/2k F ?f vM? {u) 



r(e F ,T) 2irs F Jo 1^.2 _ q 2 e(q,T) 



In the low temperature limit, with Eq. (|28p and (|29p . we find the difference between the finite 
temperature polarizability Il(g, T) and the zero temperature polarizability Il(g, T = 0) is 
just a second order (~ 0(T 2 )) small quantity. Therefore, the scattering time can be written 
as 

' ' +0(T 2 ), (45) 



t(e f ,T) t{s f ,T = 0) 
where 

1 m ^ o 2 [l-2(q/2k F ) 2 ] 2 Vl (q) 



r(e F ,T = 0) 2<Ke F J Q ,Uk 2 F -q 2 e(q,T = Q) 2 ' 



(46) 
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It is easily seen that the q ~ 2k F singularity dominate the evaluation of the integral in 
Eq. pB]) . thus we have 

1 m {2ire 2 \ 2 1 f 1 x 2 [l~2x 2 } 2 



t(e f ,T = 0) ' 27re F \ k / [1 + q TF g {2k F )/2k F } 2 Jo Vl 



dx- 



x z 



7T TLi 1 

lj £F n[C + 2k F /q TF ] 2 > ( 7) 

where n is the electron density, q FF = rage 2 / 'k is the 2D Thomas-Fermi screening wave 
vector, C = g (2k F ) = v/5 - log[(l + y/3)/2\. 

Considering the scattering time r(e, T = 0) with energy e m, we have 



71; 



2k n 2 [l-2(q/2k) 2 } 2 v t {qf 



r(e,T = 0) r (e,T = 0) n(e,T = 0) 

where 



ri(e,T = 0) vre7 2fcF ^fe 2 - g 2 ei(g,T = 0) 2 e x (q,T = 0) 

For | (s — e F )/e F |<C 1, we can write 

1 _ 1 

r (e,T = 0) ~ t(e Fi T = 0)[l + A(e - £ F )/e F ]'' 



(48) 



r(e,T = 0) 27re7 v 7 ^ 2 - q 2 e(g,T = 0) 2 ' 

where 

ei ( g ,T = 0) = l + Vc (g)iVo5o(g). (49) 
Then we express e(q, T = 0) as 

c(g,T = 0) = e 1 (q,T = 0)[1 - "'f™ ® d(q _ 2 ^)]. (50) 
With Eq. ([50j) . we can express l/r(e,T = 0) as 

1 + r7T^, (51) 



/ 2 ^ g 2 [l-2(g/2fc) 2 ] 2 Vl (q) 2 



r (e,T = 0) 27rey ^fc 2 - g 2 e 1 (?,T = 0) 

and 

1 n, /» j y[i-2(,/aff ^WiV^h) 



(54) 



with ^4 = —e F T(e F ,T = 0)d[l/T(e F ,T = 0)]/de F . From Eq. (|53l) we use the same trick as 
that of Eq. (|47p and obtain 

1 m Vi{2k F ) 2 v c (2k F )N 

:) — TTTl — ™ (55) 



n(e,T = 0) 7re v ei(2/c F ,T = 0) ' ei(2fc F ,T = 0) 
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where 



2 « q 2 [l-2(q/2k) 2 ] 2 , 2k 2 F +q 2 / 2 ., 2 . 

dq [^rv^ 2 - 4 4 + log 



2k F 



q 2 - 4k 2 F 



q + 



q 2 - Ak 2 F 



(56) 



and for | (e — £f)/ef |<S 1, we have 

2 e - e F 



+ o(( 



6 ~~ £ ^\2\ 



With Eq, ([5T1) - (f57l) . the energy dependent conductivity at zero temperature is given as 



(57) 



2„,2 



a(e,T = 0) 



e v 



e 2 v1 



-N t(e,T = 0) 
N T{e F ,T = 0) 



1 + A 



e-e F r (e,T = 0) 



e F n(e,T = 0) 



where 



r (e,T = 0) 



Aq TF 



n(e,T = 0) 2k F + Cq TF e F 
Using the Kubo-Greenwood formula [18J 



0(e-e*O 



a{e F ,T) 



1 



efe- 



cr(e,T = 0) 



cosh 1 [{e -E F )/ 2k B T] 



(58) 



(59) 



(60) 



and substituting Eq. ([58]) and ([59]) into Eq. (f60|) with a consideration of Eq. ([IT]) , we obtain 
the analytic asymptotic behavior of BLG conductivity at low temperature as following 



<t(T « T F ) = al D [l 



41og2 T 
C+l/q T F J ' 



(61) 



where qo = q FF /2k F and a 2D = e 2 VpNoT(e F ,T = 0)/2. 

For high temperature limit, substituting Eq. ([16j) into Eq. ([9]), we have 



«^>^>=^T(0( 1+c Wi)' 



(62) 



with C h = 3.57 and cr?° = (e 2 /h)(n/ni)(g 2 /irq%). 

We show our numerical results of the conductivity and analytic asymptotic result of 
Eq. (|6ip in Fig. 4, and find that they are excellently agreement with the numerical results in 
the low temperature limit. 

It is significative to compare the conductivity temperature behaviors of SLG, BLG and 
2DEG. For SLG, the asymptotic low and high temperature behaviors of conductivity are 
given by pj)] 



a(T < T F ) = a 



I- Q 



7T" 



TV 



(63) 
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Figure 4: Numerical results of temperature-dependent conductivity of BLG in the low 
temperature region and its asymptotic form given by Eq. (|6ip . The asymptotic form fits 
well in the regime T ~ 2>[0,0.04]. Here a = e 2 v 2 F N T{e F ,T = 0)/2. 



here a = a(T = 0), C ~ o(l) and J = 0.034. For 2DEG as found in Si MOSFETs and 
GaAs heterostructures, the asymptotic low and high temperature behaviors of conductivity 
are written as [19] 



here a 2D = a{T = 0), and C x = 2q Q /{l + go), a\ D = (e 2 /h)(n/m)(g 2 /irq 2 ) where q = 
q T F/2k F . 

Now let us compare the BLG temperature dependence with the SLG and the regular 
parabolic 2D systems. First, for T <C Tp, all the three systems show metallic temperature- 
dependent behaviors, but their strengthes of temperature dependence are different. BLG 
and the parabolic 2D system both have strong linear temperature dependence while SLG 
has a weak quadratic temperature dependence. Second, for high temperature limit T S> 
Tp, BLG represents a quadratic temperature-dependent behavior which is similar to SLG, 
compared with the linear temperature dependence in the parabolic 2D system. Therefore, in 
the low temperature regime, the temperature-dependent transport of BLG is qualitatively 
similar to that of the parabolic 2D system, but as the temperature increasing, BLG is getting 
more and more similarity with SLG. The transport property of BLG as the intermediate 
between SLG and the regular 2DEG has been shown here. 




(64) 




(65) 




(66) 
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4.2 Numerical Results 



We show the numerical results of resistivities obtained from Eq. Q as a function of temper- 
ature for different r s values in Fig. 5. r s = 0.88 corresponds to substrate-mounted (k = 4) 
bilayer graphene and r s = 3.51 corresponds to suspended (k = 1) bilayer graphene. It is 
found that the numerical results of BLG resistivity show metallic behavior at low tempera- 
ture and insulating behavior at high temperature, which is the same as SLG and the regular 
2D systems. Unlike graphene the scaled temperature-dependent resistivity of BLG has a rel- 
atively strong dependence on r s which is similar to that of the regular 2D systems. Because 
of the strong q = 2kp backward scattering occurring in these two systems, this similarity 
also represents in the low temperature regime, where both of them have strong linear T 
metallic behaviors with slopes 41og 2/(C + l/q ) for BLG and 2/(1 + l/g ) for 2DEG. With 
temperature increasing, the BLG temperature-dependent behavior of resistivity is changing 
from 2DEG-like to SLG-like (falling off rapidly as ~ 1/T 2 ). 

For comparison we show the calculated temperature-dependent resistivity of SLG and 
ordinary 2D systems for different interaction parameters r s in Fig. 6 (for SLG) and Fig. 7 (for 
ordinary 2D system). These two figures come from Ref. [15] . It is intuitively seen that at low 
temperature the linear T behavior of BLG is similar to that of ordinary 2D system. However, 
the linear T regime is rather weak and narrow for BLG, which is about T ~ Tp [0,0.04] for 
r s = 0.88, while the linear T regime for ordinary 2D system is relatively strong and broad, 
which extends from zero temperature to about 0.5Tp for r s = 2.6. Therefore in BLG, this 
screening-induced linear T behavior is easily suppressed by other effects. 

Since the Wigner-Seitz radius r s , representing the strength of electron-electron interac- 
tion, is reasonably small (r s ~ 0.88 for carrier density n ~ 10 12 cm~ 2 ), as shown in Fig. 5 
the temperature dependence arising from screening is rather weak in the BLG(the resistiv- 
ity of BLG for r s = 3.51 decreases just about 14 percents from T/Tp = to T/Tp = 1 
while the increase of resistivity of 2DEG for r s = 3.7 exceeds 100 percents from T/Tp = 
to T/Tp = 1). The dimensionless temperature T/Tp is also rather small because of the 
relatively high Fermi temperature in BLG (Tp ~ A00K for n ~ 10 12 cm~ 2 ). Therefore, 
as investigated experimentally in Ref. [11] and argued in Ref. [16], the strong collisional 
broadening effects due to the very small mobilities of current BLG samples would suppress 
the weak screening-induced temperature dependence that we calculated at low temperature 
(i.e., T <C Tp) and there would be no much temperature dependence in the low temperature 
resistivity. We hope our results to be tested in future experiments. 

We show the temperature-dependent conductivity of BLG for different temperatures 
calculated as a function of carrier density in Fig. 8. It shows that the conductivity increases 
in the low density regime and decrease in the high density regime as the temperature 
increases, representing a non-monotonic behavior of the conductivity. 

In this article, we have assumed an ideal 2D BLG electron gas and ignored the distance 
d between bilayer graphene and the charged impurity located at the substrate, which would 
make the form of potential of the charged impurity become 2ne 2 e~ qd /K(q + (Jtf)- We also 
have assumed a homogeneous carrier density model, and therefore our theory is quantita- 
tively correct only in the relatively high-density regime where the spatially inhomogeneous 
effects arising from charged impurity-induced electron-hole puddles are weak. For just com- 
paring theoretically the screening-induced temperature-dependent behaviors of different 2D 
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1.5 




Figure 5: Numerical results of resistivities obtained from Eq.Q as a function of temperature 
T/Tp for different r s = 3.51,0.88,0.1 (from top to bottom). r s = 0.88(3.51) corresponds 
to substrate-mounted (suspended) BLG. As r s increases the metallic behavior becomes 
stronger. 

systems {i.e., BLG, SLG and 2DEG), we do not take the level broadening effects due to 
impurity-scattering into account. 

5 Conclusions 

In this article, we calculate the static wave vector polar izability of doped bilayer graphene 
at finite temperature under the RPA. We find that for all temperatures, the BLG static 
screening is equal to its density of states Nq at zero momentum transfer and enhanced by 
a factor of log 4 at large momentum transfer. Due to the enhanced q = 2kp backward 
scattering arising from the chirality of the BLG, a strong cusp of polar izability occurs at 
q = 2fcp for zero temperature but strongly thermal suppressed as temperature increases. 
Using a microscopic transport theory for BLG conductivity at finite temperature, we also 
obtain the asymptotic low and high temperature behaviors of conductivity for BLG, and 
find it has a linear temperature metallic behavior similar to the regular 2D system at 
low temperature and a quadratic temperature insulating behavior similar to the SLG. This 
crossover from 2DEG-like behavior to SLG- like behavior as temperature increases represents 
the unique transport properties of BLG as intermediate between the SLG and the regular 
2DEG. 
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Figure 6: Numerical results of resistivities as a function of temperature T/Tp for different 
r s = 0.88, 2.2, 0.1, 0.01 (from top to bottom). r s = 0.88 (2.2) corresponds to graphene on 
the Si02 substrate (in vacuum). Inset shows the magnified view in the low temperature 
limit T < 0.5Tp. This figure comes from Ref. |15j . 
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Figure 7: p(T)/p(0) of an ordinary 2D system for different r s values as a function of 
temperature. As r s increases the metallic behavior becomes stronger. Compering to Fig. 
5, the metallic behavior of the ordinary 2D system is much stronger than that of the BLG. 
This figure comes from Ref. [15] . 
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Figure 8: Calculated conductivity for different temperatures T = 0, 100, 200, 300A" (bottom 
to top in low density regime) as a function of density. Here we use r s = 0.88 and an impurity 
density nj = 5 x W n cm~ 2 . 

References 

[1] K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. 
V. Grigorieva, and A. A. Firsov, Science 306, 666 (2004). 

[2] M. I. Katsnelson, Eur. Phys. J. B 51, 157 (2006); E. Fradkin, Phys. Rev. B 33, 3257 
(1986); A. W. W. Ludwig, M. P. A. Fisher, R. Shankar, and G. Grinstein, Phys. 
Rev. B 50, 7526 (1994); N. Shon and T. Ando, J. Phys. Soc. Jpn. 67, 2421 (1998); J. 
Tworzydlo, B. Trauzettel, M. Titov, A. Rycerz, and C. W. J. Beenakker, Phys. Rev. 
Lett. 96,246802 (2006). 

[3] E. H. Hwang, S. Adam, and S. Das Sarma, Phys. Rev. Lett. 98, 186806 (2007). 

[4] S. Adam, E. H. Hwang, V. M. Galitski, and S. Das Sarma, Proc. Natl. Acad. Sci. USA 
104,18392 (2007). 

[5] E. Rossi, S. Adam, and S. Das Sarma, Phys. Rev. B 79, 245423 (2008). 
[6] E. H. Hwang and S. Das Sarma, Phys. Rev. B 75, 205418 (2007). 
[7] T. Ando, J. Phys. Soc. Jpn. 75, 074716 (2006). 

[8] B. Wunsch, T. Stauber, F. Sols, and F. Guinea, New J. Phys. 8, 318 (2006). 

[9] Y. Barlas et al, Phys. Rev. Lett. 98, 236601 (2007). 

[10] E. McCann and V. I. Fal'ko, Phys. Rev. Lett. 96, 086805 (2006); J. Nilsson, A. H. 
Castro Neto, N. M. R. Peres, and F. Guinea, Phys. Rev. B 73, 214418 (2006); B. 
Partoens and F. M. Peeters, Phys. Rev. B 74, 075404 (2006); M. Koshino and T. 
Ando, Phys. Rev. B 73, 245403 (2006); I. Snyman and C. W. J. Beenakker, Phys. Rev. 
B 75, 045322 (2007). 



17 



[11] S. Morozov, K. Novoselov, M. Katsnelson, F. Schedin, D. Elias, J. Jaszczak, and A. 
Geim, Phys. Rev. Lett. 100, 016602 (2008); 

[12] K. Novoselov et al, Nature Phys. 2, 177 (2006); J. Oostinga et al, Nature Mater. 7, 
151 (2008). 

[13] Y. Zhang, T-T. Tang, C. Girit, Z. Hao, M. C. Martin, A. Zettl, M. F. Crommie, Y. R. 
Shen and F. Wang, Nature 459, 820 (2009). 

[14] E. H. Hwang and S. Das Sarma Phys. Rev. Lett. 101, 156802 (2008). 

[15] E. H. Hwang and S. Das Sarma Phys. Rev. B 79, 165404 (2008). 

[16] S. Adam and S. Das Sarma Phys. Rev. B 77, 115436 (2008). 

[17] S. Das Sarma and E. H. Hwang, Phys. Rev. Lett. 83, 164 (1999). 

[18] J. M. Ziman, ElectronsandPhonons (Oxford University Press, London, 1960). 

[19] S. Das Sarma and E. H. Hwang, Phys. Rev. B 68, 195315 (2003); S. Das Sarma and 
E. H. Hwang, Phys. Rev. B 69, 195305 (2004). 



18 



